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Abstract. During inflation quantum effects from massless, minimally coupled scalars and 
gravitons can be strengthened so much that perturbation theory breaks down. To follow 
the subsequent evolution one must employ a nonperturbative resummation. Starobinskii has 
developed such a technique for simple scalar theories. I discuss recent progress in applying 
this technique to more complicated models. 



1. Introduction 

On cosmological scales the universe is well described by a homogeneous, isotropic and spatially 
flat geometry, 

ds"^ = -dt^ + {t)dx -dx . (1) 

Derivatives of the scale factor a[t) give the Hubble parameter H{t) and the deceleration 
parameter 

H{t)^- , q{t)^-— = -l---^. 2 

Inflation is defined as positive expansion {H{t) > 0) with negative deceleration {q{t) < 0). 

The homogeneity of spacetime expansion evident in ([T]) does not change the fact that 
particles have constant wave vectors k, but it does alter what these mean physically. In 
particular the energy of a particle with mass m and wave number k = \\k\\ becomes time 
dependent, 

E{t,k) = ^m^+(k/a{t)y . (3) 

This results in an interesting change in the energy-time uncertainty principle which restricts 
how long a virtual pair of such particles with zizk can exist. If the pair was created at time t, 
it can last a time At given by the integral, 

t+At 

dt'E{t', fc) ~ 1 . (4) 



Just as in flat space, particles with the smallest masses persist longest. For the fully massless 
case the wave number factors out, leaving an integral which can be recognized as the radius of 
the forward light-cone from time t to time t + At [T] , 



t+At 

dt'E{t',k) 



k 



m=0 



t+At 



dt' 



a(t') q{t')H{t')a{t') 



t+At 
t 



(5) 



For positive deceleration the upper limit dominates and the integral grows without bound as 
At increases. In this case the persistence time At is finite, although longer than in flat space. 
However, for negative deceleration (inflation) it is the lower limit that dominates, and the 
integral approaches a finite value as At goes to infinity. For example, the result for de Sitter 
(a(t) = e^'' with H constant) is. 



k 



rt+At^ 
Jt a{t') 



a{t') do Sitter Ha{t) . 
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(6) 



Therefore, any massless virtual particle which happens to emerge from the vacuum with 
k < H{t)a{t) can persist forever! 

Most massless particles possess conformal invariance. The change of variables drj = dt/a{t) 
defines a conformal time in terms of which the invariant element ([I]) is just a conformal factor 
times that of fiat space, 



ds^ 



-dt^ + a^{t)dx ■ dx = —drf + dx ■ dx 



(7) 



In the [t], x) coordinates, conformally invariant theories are locally identical to their fiat space 
cousins. The rate at which virtual particles emerge from the vacuum per unit conformal time 
must be the same constant — call it P — as in fiat space. Hence the rate of emergence per 
unit physical time is, 

dN _dNdr] _ T 

~dt ~ lh]~dt ~ a(t) ■ ^ ^ 

It follows that, although any sufficiently long wavelength, massless and conformally invariant 
particle which emerges from the vacuum can persist forever during infiation, very few will 
emerge. 

Two kinds of massless particles do not possess conformal invariance: minimally coupled 
scalars and gravitons. To see that the production of these particles is not suppressed during 
inflation note that each polarization and wave number behaves like a harmonic oscillator. 



L 



1 



1 



2 2 

muj q 



(9) 



with time dependent mass m(t) = a^(t) and frequency Lo{t) = The Heisenberg equation 
of motion can be solved in terms of mode functions u{t, k) and canonically normalized raising 
and lowering operators and a, 



^2 

q + "iHq + -^q 







q{t) = u{t,k)a + u*{t,k)a^ with [a, a^] = l 



(10) 



The mode functions ti(t, k) are quite complicated for a general scale factor a(t) [2] but they 
take a simple form for de Sitter, 



u(t, k) 



H 

7W 
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ik 
Ha{t). 



exp 



ik ' 
Ha{t). 



(11) 



The (co-moving) energy operator for this system is, 



E{t) = lm{t)q\t) + ^m{t)u;'{t)q\t) . (12) 

Owing to the time dependent mass and frequency, there are no stationary states for this system. 
At any given time the minimum eigenstate of E{t) has energy ^u>{t), but which state this is 
changes for each value of time. The state \^}) which is annihilated by a has minimum energy 
in the distant past. The expectation value of the energy operator in this state is. 



(n E{t) n) = la^{t)\u{t, k)\' + la{t)kMt, =^ + ^. (13) 

\ / 2 2 de Sitter la 4fc 

If one thinks of each particle having energy k/a(t), it follows that the number of particles with 
any polarization and wave number k grows as the square of the inflationary scale factor. 

Quantum field theoretic effects are driven by essentially classical physics operating in 
response to the source of virtual particles implied by quantization. On the basis of ()14p 
one might expect inflation to dramatically enhance quantum effects from MMC scalars and 
gravitons, and explicit studies over a quarter century have confirmed this. The oldest results 
are of course the cosmological perturbations induced by scalar inflatons [3] and by gravitons |4j. 
More recently it was shown that the one loop vacuum polarization induced by a charged MMC 
scalar in de Sitter background causes super-horizon photons to behave like massive particles 
in some ways [Sj El [7j . Another recent result is that the one loop fermion self-energy induced 
by a MMC Yukawa scalar in de Sitter background reflects the generation of a nonzero fermion 
mass [SI E] • In the next three sections it will be explicitly shown how these one loop results 
generalize to all orders. 

2. Infrared Logarithms 

The expectation values of familiar operators typically show enhanced quantum effects in the 
form of infrared logarithms. A simple example is provided by the stress tensor of a massless, 
minimally coupled scalar with a quartic self-interaction, 

1„ „ „., , A 



k H^a 



C = —-d^ipdijipg^'^yj—g — -^'P^V~9 + counterterms . (15) 

When the expectation value of the stress tensor of this theory is computed in de Sitter 
background {a{t) = e^*) and renormalized so as to make quantum effects vanish at t = 0, 
the results for the quantum- induced energy density and pressure are [lOl [H] , 



Pit) = j^,{-ln\a)}+0{X'), (16) 



(2^) 

P^'^ = (^{-^^^'^"^-1^^"(")} + ^^^')- (^^^ 

Infrared logarithms are the factors of ln(a) = Ht. They arise from the fact that inflationary 
particle production drives the free scalar field strength away from zero [12 ^ 113 ^ [T^. 



2 



(/?^(x) 17) = ln(a) + Divergent Constant . (18) 



This increases the vacuum energy contributed by the quartic potential, and the result is evident 
in ^M3). 

Infrared logarithms arise in the one particle irreducible (IPI) functions of this theory [15j . 
They occur as well in massless, minimally coupled scalar quantum electrodynamics (SQED) 
O El [7] and in massless Yukawa theory [HJ [9] . The IPI functions of pure gravity on de Sitter 
background show infrared logarithms [161 I17j . It seems inevitable that infrared logarithms 
contaminate loop corrections to the power spectrum of cosmological perturbations [181 [T9] and 
similar fixed-momentum correlators ^20j. And infrared logarithms have been discovered in the 
IPI functions and quantum-corrected mode functions of Einsetin -|- Dirac |21[ I22j . 

Infrared logarithms are fascinating because they introduce a secular element into the usual, 
static expansion in the loop counting parameter. No matter how small the coupling constant A 
is in (|16til7p . the continued growth of the inflationary scale factor must eventually overwhelm 
it. When this happens, perturbation theory breaks down. For example, the general form of 
the induced energy density p6|) is, 

oo 

p(t) = H''J2^'''{Co^^''~Ha) + Ciln^'-'ia) + ... + C',,_,ln\a)] . (19) 
e=2 

The CQ[Aln^(a)]^~^ terms are the leading logarithms at £ loop order; the remaining terms are 
subdominant logarithms. Assuming that the numerical coefficients C| are of order one, we 
see that the leading infrared logarithms all become order one at ln(a) ~ l/^/X. At this time 
the highest subdominant logarithm terms are still perturbatively small (~ \/A), so it seems 
reasonable to attempt to follow the nonperturbative evolution by resuming the series of leading 
infrared logarithms, 

Plead = ^^5Z'^o(^ • (20) 

£=0 

This is known as the leading logarithm approximation. 



3. Starobinskii's Formalism for Simple Scalar Models 

Starobinskii has long maintained that his stochastic field equations reproduce the leading 
logarithm approximation [23]. With Yokoyama [21] he exploited this conjecture to explicitly 
solve for the nonperturbative, late time limit of any model of the form , 

C = -^df,^d,^g^''^-V{^)^ , (21) 

assuming only that the potential V{(p) is bounded below. When the potential is unbounded 
below the conjecture still gives the leading infrared logarithms at each order, however, the 
theory fails to approach a static limit. 

An all orders derivation has recently been given of Starobinskii's formalism [25^ I26j . The 
first step is to rewrite the operator field equations, 

in Yang-Feldman form |27j . 

ip{t, x) = ifoit, x)- dt'a'^ J d^x'G{x; x')V' ((^(x')) . (23) 



(22) 



The free field expansion and retarded Green's function are, 



G{x;x') 



(27r)3 
ie{At) 



u{t, A;)e*^-^a(fc) + u*{t, A;)e-^^ V(fe)} 
(2^ 



Ak-Ax 



^u{t,k)u*{t',k) - u*{t,k)u{t',k)^ . 



(24) 
(25) 



The mode function u{t,k) was given in expression (llip . The canonicahy normahzed creation 
and annihilation operators are a^{k) and a{k). Iterating the Yang-Feldman equation gives the 
usual perturbative expansion of the interaction picture field, expressed in terms of a field which 
is free at t = 0. 

Now consider taking the expectation value of some operator constructed from ip, and hence 
from the free field (fQ . To reach leading logarithm order requires that every free field contributes 
to an infrared logarithm. The full result from the pairing of two free fields is, 



{t,x)ipQ{t' ,x) Qq^ = 1 dkk'^ u{t,k)u*{t' ,k) 



, sin(A;Ax) 



kAx 



(25) 



At high k the mode functions and the sin(fcAx) oscillate, which makes the integral converge. 
As the name suggests, it is the low k end of the integration which is responsible for infrared 
logarithms. In this regime sm(kAx)/kAx ~ 1 and only the first term in the long wavelength 
expansion of the mode functions matters. 



u{t, k) 



H 



/2P 



(26) 



This observation has two important consequences: 

• The leading logarithm result will not be changed if the free field mode sum is cut off at 

k = Ha{t); and 

• The leading logarithm result will not be changed if the mode function is replaced by its 
infrared limit. 

Together, these simplifications convert the free quantum field to a commuting combination of 
creation and annihilation operators. 



I ^9{Ha - k)-^{e^'Mk) + e-^^-V(^)} , (27) 



Because the Green's function involves a commutator of mode functions, its infrared truncation 
requires third order terms from ()26|) . 



G{x; x) 



ie{t - t') 

e{t-t') 

3H 



A;3 \3\HaJ 3\Ha'J i ' 



(27r)3 



a'3 



5^{x-x') 



The infrared truncated Yang-Feldman equation is accordingly. 



^t,x) = Mt,x)- ^ I dt' 



1- (- 

a 



V'(^{t',x) 



(28) 
(29) 

(30) 



Although ip{t,x) and ^{t,x) are vastly different operators, expectation values involving them 
agree at leading logarithm order. 

If we assume the scalars in V' (^<p{t', x)^ grow like In(a') — which is certainly true whenever 

VEV's are taken — then one sees that the t' integration in ()30p can produce an additional 
infrared logarithm for the first of the square-bracketed terms. However, the rapid growth 
of the term (a'/a)^ weights the integral overwhelmingly at its upper limit and precludes the 
development of an additional infrared logarithm. We can therefore ignore this term and simplify 
to the equation, 

1 



Taking the time derivative gives Starobinskii's Langevin equation [24], 

Mt,x) = ^o{t,x) - (<!>{t,x] 



(31) 



(32) 



Starobinskii's stochastic noise term is the time derivative of the infrared truncated free field, 



i>o(t,x) 



(27r)3 ^ '^2k^ 
A simple calculation reveals that it behaves like white noise 



n 



(33) 



(34) 



Langevin equations of the form (j32p have been much studied [28j. Expectation values of 
functionals of the stochastic field can be computed in terms of a probability density Q{t, (j)) as 
follows, 



F[4>{t,x)] 



d4>Q{t,4>)F{(t>) 



(35) 



The probability density satisfies a Fokker-Planck equation whose first term is given by the 
interaction in (j32p and whose second term is fixed by the normalization of the white noise p4p : 



1 a' 

^ 25^ 



|2 r t73 



(36) 



To recover the nonperturbative late time solution of Starobinskii and Yokoyama [23] one makes 
the ansatz, 

lim £>(t,0) = £»oo(0) , (37) 

t— »oo 

because the scalar force should eventually balance the tendency of inflationary particle 
production to force the scalar up its potential. This ansatz results in a first order equation, 



dQc 



Qoo{4>) 



8^ 
3i^4 



y'(0)# . 



(38) 



The solution is straightforward. 



Qoo{4>) = Nexp 



Svr^ 



(39) 



Figure 1. Two loop contribution to {Q\Fn^{x)Fp(^{x)\Q) . 




Figure 2. Effective {(p*^p)'^ coupling in SQED. 
4. More General Scalar Models on de Sitter 

A field which can generate infrared logarithms is called active. Scalar potential models of the 
form (j2ip possess only active fields. However, more general theories can possess passive fields 
which are not themselves capable of engendering an infrared logarithm. A example of such a 
model is scalar quantum electrodynamics (SQED), 

(40) 

In this model the charged scalar is active whereas the photon is passive. 

Although passive fields cannot cause infrared logarithms, they can propagate their effects. 
That is, an expectation value of passive fields can acquire an infrared logarithm from a loop 
correction involving an active field. For example, the diagram in Fig. 1 gives a contribution 
to {fi\Ffj_i^{x)Fp(j{x)\Q) which acquires an infrared logarithm through the scalar loop at the 
bottom. Passive fields can also induce interactions between active fields. For example, the 
photon loop in Fig. 2 induces an effective (95*93)^ interaction in SQED. 

In generalizing Starobinskii's technique to theories which include passive fields, it is crucial 
to realize that the ultraviolet parts of passive fields fields contribute on an equal footing with the 
infrared parts to the processes of propagating infrared logarithms and mediating interactions 
between active fields. So one cannot infrared truncate the passive fields. Instead the correct 
procedure is: 

• Integrate out the passive fields and renormalize the resulting effective action; then 

• Infrared truncate and stochastically simplify the purely active field effective action. 

One might suspect that the second step is not possible owing to the nonlocality of the effective 
action. However, first note that VEV's have no spatial dependence at leading logarithm order, 
so one can evaluate all fields at the same space point. It is still necessary to confront the 
prospect of fields buried inside different temporal integrals. However, precisely because they 
derive from passive fields, these temporal integrations always contain positive powers of the 



scale factor whose rapid time dependence weights the integral overwhelmingly at its upper limit 
and totally dominates the logarithms which might derive from the active fields. 

As an example, consider the retarded Green's function of a conformally coupled scalar, 



Gef(x,x ) - ^e{At) . (41) 



Acting this Green's function on In(a') gives. 



dt'a'^ I d^x'Ga{x]x') X In(a') = — f dt'\- - {-\^] In(a') , (42) 

H Jq V a \ a / \ 



2i?2 

(To appreciate the distinction between the Green's functions of passive and active fields, 
contrast (|42|) with (j3U|) .) On the other hand, the result of simply multiplying ln(a) by the 
Green's function acted upon unity is, 

ln(a) X dt'a'^ j d^x'Grf (x; x') = ^{l - ^ + ^} • (44) 

These expressions agree at leading logarithm order, so one may as well locate all the scalar 
fields at the same point and evaluate the inverse differential operators on unity. But that is 
just the same thing as computing the effective potential! Hence the hopelessly complicated 
"effective action" degenerates, in the leading log approximation, to a very tractable "effective 
potential," and the resulting local theory assumes the form (12 ip already solved by Starobinskii 
[231 [2i]. 

This program has been carried out for SQED in collaboration with Nikolaos Tsamis and 
Tomislav Prokopec [29]. The various results are best expressed in terms of the following 
quantity, 

z = . (45) 

If one renormalizes to make the quadratic and quartic terms of the effective potential vanish, 
the result is, 

^«ff = ^|(-l + 27)^+ (-^+7)^' + ^^ dx(l+x)[v42-Az.o) + V'(l + Aj^o)]| • (46) 

Here the function IS.vq{x) is, 

Ai/o(x) = i - ivi-8x = 2x + 4x2 + 8x^ + 0{x^) , (47) 
and the PolyGamma function is. 



^{l + z) ^ -ln(r(l + z)) = -^ + ^(-l)X(n)z"-i . (48) 

n=2 



The effective potential of SQED in de Sitter background does not seem to have been computed 
previously. However, the H = limit of expression ()46p agrees with equation (4.5) of Coleman 
and Weinberg [30] . 



To evaluate the VEV of any operator, one first integrates out the passive fields and 
stochastically simplifies the resulting, purely active field functional. One then computes the 
VEV using Starobinskii's formalism. Because the ultraviolet contributes for passive fields, the 
VEV's of some operators are ultraviolet divergent even at leading logarithm order. The scalar 
functional resulting from the Z)-dimensionally regulated field strength bilinear is, 



ffD r(f-i)r(2-f)r(^+i/)r(^ 



(4vr: 



r(#+i) 



Here the scalar-dependent parameter 2^(2;) is. 



L>-3\2 /D-3\ , , , 

^) -2.-(^)-AK^). 

The analogous result for the scalar kinetic term is, 

HD (r{§-l)T{2-§)T{^ + u)T{^-i^) r{D-l) 



4 (47r) 



r(f+i) T{^^+u)T{i-u) r(f+i) 



(49) 



(50) 



(51) 



The leading logarithm result for the stress tensor is finite and takes the form —g^i,Vs{z) where. 



16^ 



-2+27)z+ [-^ + 2^]z^ + z{l+z) V'(l-2Az^o) +V'(l + Ai^o) 



(52) 



The stochastic prediction for (|51|) has been checked by an explicit two loop computation |31| . 
The predictions for (|51|) and (j52|) are still being checked. 

In collaboration with Shun-Pei Miao, Starobinskii's formalism has also been applied to a 
massless, minimally coupled scalar which is Yukawa-coupled to a massless fermion [32J. The 
Lagrangian of this model is. 



1 



1 



g + i■^l)e\-i^Vp%l)^f^ - f^p'^pil^f^ , (53) 



where the spinor covariant derivative is. 



9fi + —A.pcdJ'^'^ 



(54) 



As with SQED, the conformal and quartic counterterms can be chosen to make the quadratic 
and quartic terms in the renormalized effective potential vanish. The final result is. 



y> (-1)" 
47r2 ^ n+1 

n=2 



C(2n-1) -C(2n+1^ 



(55) 



H J 



+ 2 1 dx{x+x^) ^j{l+ix)+%lj{l-ix) ). (56) 



This agrees with the classic result of Candelas and Raine [33l [Ml [35] , and of course the H = 
limit agrees with equation (6.10) of Coleman and Weinberg [30l[36]. Note that this potential 
is unbounded below, which establishes that infrared logarithms need not always sum up to 
approach a static limit at late times. 

As was the case for SQED, the leading logarithm result for the Yukawa stress tensor takes 
the form, —g^^Vs{if), where 



^-7+C(3) 



SlKif)^+Kif)*-E(-r[c(2„-i)-c(2„.i) 



87r2| 



n=2 




The leading logarithm result for the dimensionally regulated coincident vertex function is, 

2 



(p{x)ij){x)ijj{x) 



H 



D-2 



(47r)f 



■ 2 ^ H ' 



r(i + ^ 



H 



Expression (|59|) was checked at order / against an explicit two loop computation 



(59) 



5. Other Geometries 

Applying Starobinskii's formalism to more general theories on de Sitter background is essential 
to resolve the issue of what happens when infrared logarithms in these models become 
nonperturbatively strong. However, this does not suffice for extrapolating the enhanced 
quantum effects of these models to post-inflationary cosmology. For that it is necessary to 
understand how the effects manifest for a general scale factor a{t). Two issues are of special 
importance: 

• What becomes of the dependence upon the inflationary Hubble parameter H which arises 
from integrating out passive fields? 

• What becomes of the infrared logarithms generated by active fields? 

Because the passive fields give local effects, it seems plausible that factors of H generalize 
to local curvature scalars. The unique curvature scalar of dimension two is the Ricci scalar. 
Evaluating it respectively for a general scale factor and for de Sitter gives. 



12^2 (t) + m{t) — > 12H^ (de Sitter) . 



(60) 



This suggests that factors of in de Sitter generalize to R/12. For example, the 
generalizations of the effective potentials of SQED PB]) and Yukawa would take the form, 



SQED 



2 * 

e ip (f 

/V 



(61) 
(62) 



It is presumably this additional dependence upon the metric which is responsible for the fact 
that the potentials of the leading logarithm stress tensors — expressions (I52p and ()57p — do 
not agree with the effective potentials (I46p and (I56p . Note also that the putative generalization 
produces a curious sort of R'^f{ip'^/R) model which may itself have cosmological significance. 

Even in de Sitter, infrared logarithms measure the time since the onset of inflation, so they 
cannot generalize to a local invariant. One way of inferring the nonlocal invariant to which they 
generalize is to consider the expectation value of the square of a free, massless, and minimally 
coupled scalar for an arbitrary scale factor. Because the result is ultraviolet divergent, it must 
be regulated, although the generalized infrared logarithm should be finite. If the general a(t) 
mode function is denoted u[a](t, k), the result can be written as follows, 



Of 



ipl{t,x] 



Of 



d^-^k 



(2vr) 



D~l 



u[a]{t, k) 



(4vr)^r(^) Jh 



dkk 



D-2 



u[a]{t, k) 



(63) 

Explicit expressions for u[a]{t,k) exist [2l[37j. Although these expressions are very complicated, 
the leading time dependence of (fH5|) can perhaps be extracted. 



6. Discussion 

There is no question that infrared logarithms arise in explicit perturbative computations 
in many inflationary quantum field theories which involve MMC scalars and/or gravitons 
[SlElElElElIinillllliailllllTlIIElEIlISn^ These infrared logarithms counteract the 

small coupling constants which would otherwise suppress quantum loop effects. Over a long 
period of inflation they become so large that weak field perturbation theory breaks down. 
When the series of leading infrared logarithms is summed using Starobinskii's formalism, the 
result is that MMC scalars always reach nonperturbatively large field strengths. In the case 
of Yukawa theory the scalar grows without bound and comes to dominate late time cosmology 
[32j . It is not yet known what the nonperturbative outcome is for quantum gravity. 

These enhanced quantum effects are not restricted to the MMC scalars and gravitons which 
cause them. Other fields seem to experience the following effects: 

• Scalars with nonderivative couplings induce growing masses [SI El [TJ [U El [IS] ; and 

• Gravitons and scalars with derivative couplings seem to induce changing field strengths 

[lElES]. 

The now-established fact of significant, nonperturbative quantum effects from MMC scalars, 
and the prospect for them from gravitons, change the inflationary paradigm. Four processes 
occur which are absent in classical inflation: 

• Certain fleld strengths reach nonperturbatively large values (~ /coupling const); 

• Inflation-induced masses engender a signiflcant amount of vacuum energy which is positive 
for vector bosons and negative for fermions; 

• The effective action of gravity suffers modifications of the form A£ = E? f {ip^ / R); and 

• A potentially significant amount of negative vacuum energy derives from the self- 
gravitation of inflationary gravitons. 

Each of these processes can produce observable signals: 

• Evolving scalars can induce resonant particle production [55] : 



• Vaccum energy from (Hubble scale) particle masses can serve as a reservoir for the creation 
of a vast population of cosmological wavelength vectors [39l HOj HI] and fermions when the 
mass dissipates after second horizon crossing; 

• Modifications of gravity can change cosmology [321 IS] ; and 

• It is conceivable that inflation was caused by a bare, GUT-scale cosmological constant 
which is being screened by infrared quantum gravitational effects that became strong 
during primordial inflation [44j . 
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